Band-like transport and trapping in TMTSF Single Crystal Transistors 
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We perform a combined experimental and theoretical study of tetramethyltetraselenafulvalene 
(TMTSF) single-crystal field-effect transistors, whose electrical characteristics exhibit clear signa- 
tures of the intrinsic transport properties of the material. We introduce a simple, well-defined model 
based on physical parameters and we successfully reproduce quantitatively the device properties as 
a function of temperature and carrier density. The analysis allows its internal consistency to be 
checked, and enables the reliable extraction of the density and characteristic energy of shallow and 
deep traps in the material. Our findings indicate that shallow traps originate from electrostatic 
potential fluctuations generated by charges fixed in the deep traps. 
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Despite impressive progress in the use of organic semi- 
conductors for the realization of practical circuits and de- 
vices, our fundamental understanding of these materials 
remains limited. In this context, the recent observation 
of an anisotropic carrier mobility increasing with lower- 
ing temperature in organic single-crystal field-effect tran- 
sistors (FETs)pj establishes a breakthrough, because it 
opens the possibility to investigate intrinsic transport 
properties of molecular materials at finite and tunable 
carrier densityQ, Surprisingly, a systematic, in-depth 
analysis of the intrinsic transport phenomena observed 
in single-crystal FETs has not been attempted yet. In 
large part, this is due to the lack of a well-defined model 
enabling a quantitative study to be performed. In fact, 
past work on the intrinsic transport properties of organic 
semiconductors -based on time-of-flight measurements- 
has relied on the so-called trap-and-release model 0, Hj], 
which, despite being helpful 6], cannot be applied to the 
systematic analysis of systems at finite carrier density. 
This poses a clear problem, because a quantitative analy- 
sis of the properties of organic semiconductors is essential 
to characterize devices and materials unambiguously, to 
perform comparative studies, and to reach a true micro- 
scopic understanding, which is essential for the efficient 
long-term development of practical applications. 

We address this issue through a combined experimen- 
tal and theoretical study of transport through tetram- 
ethyltetraselenafulvalene (TMTSF, Ci Hi 2 Se4, see Fig. 
la) single-crystal FETs. We first show that these de- 
vices exhibit the characteristic signature of the intrinsic 
transport properties. We then introduce a simple model 
based on physically transparent quantities and assump- 
tions, with which we analyze in detail the temperature 
and density dependence of the mobility, and the temper- 
ature dependence of the threshold voltage. Our analysis 
reproduces the observations systematically and quantita- 
tively, it gives consistent results on different devices, it 
enables the reliable determination of the concentration 
and characteristic energy of traps, it allows the check of 
its internal consistency to be performed, and it points 
to the physical origin of the shallow traps. These results 



clearly indicate that the model represents an appropriate 
framework for the analysis through high-quality organic 
FETs. 

So far, rubrene crystals have been the only ones 
to unambiguously show the expected signatures of in- 
trinsic transport at finite carrier density in a FET 
configuration[l|. We have explored crystals of different 
molecules, and here we discuss TMTSF transistors 0, re- 
alized by laminating a single-crystal grown from vapor 
phase transport onto a polydimethylsiloxane (PDMS) 
support as shown in Fig. lb,c (the fabrication proce- 
dure is identical to that described in Ref. [8J)- Immedi- 
ately after fabrication, the devices were transferred into 
the vacuum chamber (p w 2 10~ 7 mbar) of a flow cryo- 
stat where all transport measurements were done. Most 
measurements were made using a HP 4156A parameter 
analyzer, in a four-terminal configuration (Fig. 1 (c)), to 
eliminate the effects of the contact resistance. We discuss 
data taken on three different FETs, representative of the 
approximately 10 devices studied. 
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FIG. 1: (a) Chemical structure of TMTSF molecule, (b) 
Schematic cross section of our PDMS-based TMTSF transis- 
tors, (c) Optical microscope image of a TMTSF PDMS FET 
(the seperation between source and drain is 300 fj,m). 

Fig. 2(a) shows the source-drain current (Isd) mea- 
sured at different temperatures, as a function of gate volt- 
age ( Vg ) , and the inset Is d measured on the same device 
as a function of Vsd at T=300 K, for different values of 
Vq- The measurements exhibit the characteristic behav- 
ior of high-quality devices (e.g., absence of hysteresis, 
stability and reproducibility). The charge carrier mo- 
bility is extracted from the linear part of Isd-Vg char- 
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acteristics, using the relation \i — ^tTF cWb (with L 
separation between voltage probes, W crystal width, d 
gate capacitance per unit area, and V voltage difference 
measured between the two voltage probes) . At room tem- 
perature, fj, = 4 cm 2 /Vs reproducibly (the same value is 
obtained from the Isd — Vsd characteristics in the satu- 
ration regime) . From the Isd ~ Vg curves we extract the 
full temperature dependence of the mobility, as well as 
the threshold voltage Vt (T) (by extrapolating the curves 
to zero current). 
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FIG. 2: Source-drain current, measured as a function of gate 
voltage at different temperatures (Vsd=1 V). The inset shows 
the source-drain current as a function of source-bias voltage, 
for different values of Vg (from -IV to -19V, in -2 V steps). 

Fig. 3a and b show the mobility and threshold voltage 
as a function of temperature measured on three differ- 
ent samples. The mobility tends to increase, exhibiting 
an identical temperature dependence down to a sample- 
dependent temperature value, after which /z decreases 
with further lowering T. The temperature interval where 
fi(T) is "metallic- like" reflects the intrinsic material prop- 
erties, and the different extents of this interval in the 
three samples indicate that different amount of disorder 
is present. Interestingly, the decrease of the mobility at 
low T is very slow. In our best sample, ji reaches a value 
above 6 cm 2 /Vs at T ~ 160 K, and it is still close to 2 
cm 2 /Vs at 50 K. The behavior of the threshold voltage 
is qualitatively similar in different devices -exhibiting a 
linear decrease with lowering T at first, and a tendency 
to saturate at lower T. The magnitude of Vt depends 
on the sample, and it is systematically smaller in higher 
mobility devices (i.e., the experiments indicate that the 
temperature dependence of fi correlates with the values 
of Vt). 

To analyze quantitatively the data we introduce a 
model based on the assumption that transport of holes in 
TMTSF occurs in a band (with density of states No for 
energy E > 0), and that disorder- induced trap states are 
present (for E < 0)(see also 9]). These consist of shal- 
low traps (with density of states Ns{E)) and deep traps 
located much deeper in energy (described by a density of 



FIG. 3: The symbols represent the temperature dependent 
mobility (a) and threshold voltage (b) measured on three dif- 
ferent devices. The solid lines represents the result of our 
model, using the parameters listed in Table 1; the red and 
black lines are obtained using a constant and a gaussian den- 
sity of shallow traps, respectively. 



states Nd(E)). The mobility of carriers occupying states 
in the band is the intrinsic mobility in the material, and, 
following existing theories, is assumed to depend on tem- 
perature as hq(T) = aT- 2 0, [Til [l"2T | . Carriers occupy- 
ing trap states are assumed not to carry current. 

In practice, we fix the shape of the density of states and 
calculate the position of the Fermi level Ep (numerically) 
as a function of the total density of accumulated charge 
n (= CiVc) through the relation: 
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that defines Ep. Having determined Ep(n,T), it is 
straightforward to obtain the source-drain current Isd in 
terms of the density of charge carriers occupying states 
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From Eq. 2 we extract all quantities that are measured in 
the experiments, by applying to the calculated Isd — Vsd 
curves the same procedures used for the measured ones. 
For instance, the mobility \i measured from the FET 
characteristics is simply given by fj, = (J-q{T) d" g^'^ 
(which corresponds to fj, a d gy® ) ■ For our calculations 
we use, for the shallow traps, a constant density of states 
(N S (E) = N c = const, if E c < E < 0, and N S (E) = 
otherwise), and a gaussian density of states (Ns(E) = 
N g e~ l y E l 2E 9) ) 5 and compare the results obtained impos- 
ing that the total number of shallow traps /_ Ns(E)dE 
is the same for the two distributions. For the deep traps 
we use a "square" distribution (Nd(E) = Nd = con- 
stant if Ed - SE D /2 < E < E D + SE D /2 and zero 
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otherwise) to analyze the role of the trap depth (Ed) 
and of the distribution width (AEjj). The density of 
states in the band at the surface of the TMTSF crystal - 
N = 10 15 cm -2 eV _1 - is estimated by taking one state 
per molecule distributed in energy over the bandwidth 
of the valence band (~ 0.5 eV [13|). The parameter a 
(in the expression Hq(T) — aT~ 2 ) represents an intrinsic 
property of TMTSF crystals, and it is taken to be the 
same for all samples. 

A strategy similar to the one adopted here -i.e., "fix- 
ing" a density of states to analyze transport through 
FETs- is used to describe transport in FETs based on 
amorphous materials In that case, however, trans- 
port is due to carriers occupying localized states, whose 
mobility is strongly energy dependent, and material in- 
homogeneity requires percolation effects to be taken into 
account [lij . This significantly increases the complexity 
of the analysis, as well as the unambiguous extraction of 
parameters in the model. These problems are absent in 
high-quality single-crystal FETs. Only at low temper- 
atures, such that most carriers occupy trap states, the 
model starts to deviate from the experiments, mainly be- 
cause our assumption that carriers in the shallow traps do 
not carry any current is too restrictive. For this reason, 
we limit our analysis to a temperature range where a sig- 
nificant fraction of carriers populate states in the band. 
In this temperature range (between approximately 100 
and 300 K for the best device) we find that the behavior 
of the temperature dependence of the mobility /i(T) mea- 
sured at high carrier density curve depends only on the 
properties of the shallow traps (Fig. 3a), whereas the 
behavior of the Vr(T) is determined only by the deep 
traps. This considerably simplifies the analysis. 




FIG. 4: In panel (a) and (b) the symbols show the mobility 
extracted from FET measurements as a function of carrier 
density at different temperatures, on devices 1 and 2. The 
solid lines are calculated using the parameter values extracted 
from the fits shown in Fig. 3 (no other parameters need to be 
introduced). For device 3, a similar agreement between data 
and model is obtained. 

Best fits to the temperature dependence of the mobil- 
ity are shown in Fig. 3(a), using both the gaussian and 
the constant distributions of shallow traps. With both 
distributions the agreement is satisfactory for the tem- 
perature range discussed above (ljj. Table 1 summarizes 



the parameters extracted from the fits. Note that the 
energies E c and E g satisfy the relation E c = 3E g to a 
10% precision (or better) for all samples, which can be 
easily understood since the integral of the gaussian trap 
distribution between E = 3E g and E = gives 99% of 
the shallow trap states. Since different distributions lead 
to nearly identical temperature dependencies of the mo- 
bility, the precise energy distribution of the shallow traps 
cannot be extracted from the analysis. On the contrary, 
the total density of shallow traps and of their character- 
istic energy (~ E g ) are not sensitive to the details of the 
model, and the values extracted have therefore physical 
meaning. 

The microscopic analysis of the threshold voltage is 
meaningful, because in our FETs the single-crystals are 
suspended on top of the gate, which eliminates extrinsic 
effects due to phenomena taking place in the gate insu- 
lator. We find that the curve Vt(T) is sensitive to both 
the energy and the width of the deep trap distribution. 
Specifically, very deep traps (whose energy is larger than 
roughly 0.5 eV) contribute to the threshold voltage, but 
do not cause a shift in the temperature range considered. 
In this same range, a sizable contribution to the threshold 
voltage shift originates from traps with energy between 
100 and 400 meV, approximately, with the largest con- 
tribution coming from traps between 150-200 meV. This 
fact -that deep traps at different energies affect differ- 
ently Vt and its shift with temperature- makes it pos- 
sible to extract information on the energy distribution 
of deep traps from the experimental Vt{T) curves. With 
the parameter values listed in Table 1, the model leads to 
an excellent quantitative agreement with the data (again, 
below approximately 100 K, where the majority of car- 
riers occupies shallow traps, the results depend on the 
details of the shallow trap distribution). Note that the 
center of the deep trap distributions is approximately the 
same in the three devices. This is expected if deep traps 
originate from a well-defined chemical impurity hosted 
in the crystal. Indeed, the values of Ed and of SEd ex- 
tracted from our analysis are in the same range energies 
and widths of other deep traps recently found in different 
molecular materials (l7). 

With the model reproducing successfully the exper- 
imental observations, it is important to critically ask 
whether the quantitative agreement is not simply due 
to the introduction of a sufficient number of free parame- 
ters. To show that this is not the case, we have analyzed 
the full dependence of the mobility on carrier density, 
at different temperature values. This comparison pro- 
vides a stringent test of the validity of the model and of 
its consistency, because all parameters are fixed by the 
analysis discussed above, and no other free parameters 
can be introduced. Fig. 4a and b show that the curves 
calculated with the model are in excellent quantitative 
agreement with the data, as long as the temperature is 
such that enough charge carriers are present in the band. 
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S ample 1 


Sample2 


Sample3 


N c (cm- 2 eV" 1 ) 
Ec(meV) 


1.97*10 14 
38.5 


8.5*10 13 
54.1 


9.3*10 12 
111 


N g (cm- 2 eV" 1 ) 
E 9 (meV) 


3.04*10 14 
14.1 


1.36*10 14 
19.1 


1.42*10 13 
40.8 


J° oo Ns(E)dE (cm- 2 ) 


7.6*10 12 


4.6*10 12 


1.03*10 12 


Ni,(cm- 2 eV _1 ) 

E D (meV) 

SE D {meV) 


6*10 10 

225 

95 


1.08*10" 

235 

95 


1.08*10" 

270 

130 


A$(meV) 


10.3 


13.7 


16.1 



TABLE I: Values of quantities discussed in the text. The total 
number of shallow traps N,s (E)dE is fixed to be the same 
for both the constant and the gaussian energy distributions. 

This result confirms the internal consistency of our anal- 
ysis, and illustrates that the model does describe well the 
interplay between band-like transport and trapping as a 
function of carrier density. 

Having established the internal consistency of our anal- 
ysis, we address the origin of the observed correlation 
between the temperature dependence of mobility and of 
threshold voltage. Since the temperature dependence of 
the mobility is determined only by the shallow traps and 
that of the threshold voltage only by the deep traps, this 
correlation implies that a physical mechanism exists, that 
links the origin of the two kinds of traps. We propose 
that this mechanism is due to simple electrostatics. In- 
dividual charge carriers frozen in the deep traps gener- 
ate random electrostatic potential fluctuations A<I>, and 
the resulting electrostatic potential landscape causes re- 
gions of lower energy in which a large density of charge 
carriers can be trapped. The characteristic energy scale 
associated to the potential fluctuations can be estimated 
as A<£> = -r-= — -zr^-, where < r > is the average distance 

47reeo <r> ' G 

between deep traps and e is the average of the relative di- 
electric constant of the TMTSF crystal and of vacuum. 
By estimating the magnitude of < r > from the total 
density of deep traps extracted from the fit of Vt(T), we 
obtain the values of A<!> reported in Table 1 , whose mag- 
nitude compares well with the characteristic energy scale 
Eg of the shallow traps. This scenario explains why the 
density of states of shallow traps at E ~ is comparable 
to the density of states in the band itself. In fact, shallow 
traps consist of states in the band whose energy is low- 
ered by the local, slowly varying electrostatic potential. 
It also explains why carriers in shallow traps contribute 



significantly to the conductivity (as we noticed earlier), 
since these carriers are not tightly bound to a defect and 
can move over rather large distances. 

In conclusion, we have developed a well-defined frame- 
work to analyze the interplay between intrinsic band-like 
transport and trapping in organic single-crystal transis- 
tors and used it to describe the behavior of newly devel- 
oped TMTSF FETs. The model permits to understand 
in detail the physical origin of the different aspects of the 
device characteristics, and to extract microscopic param- 
eters reliably. This combination of experimental investi- 
gation on high-quality devices and systematic quantita- 
tive analysis is important to reach a true understanding 
of molecular semiconductors at finite density of charge 
carriers. 
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